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Abstract 

We consider in a pedagogical fashion alterations to Newtonian gravity due to the 
postulate that all energy corresponds to active gravitational mass when applied 
to the self-energy of the gravitational field. We show why a simple addition of 
4j- times the gravitational field energy to the matter density in Newton's field 
equation is inconsistent. A consistent prescription is shown and discussed. The 
connection to general relativity is pointed out. 



1. Introduction 

The issue addressed in this letter arises if one wishes to model the self-coupling of 
the gravitational field within Newtonian gravity. Simple non-linear alterations of 
Newton's field equation are often employed as simplified models for general relativ- 
ity. The purpose of this letter is to show how this can be done and to point out 
certain flaws in the usually accepted prescription, as for example given in [1] [2] . 

To be more precise, we recall that the Newtonian gravitational field, (p, and 
the density of (ponderable) matter, p, obey 

Aip = 4nGp (1.1) 

where G is Newton's gravitational constant. The force per unit volume is given by 

/ = -pV<p. (1-2) 
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Together these equations imply that in order to build up a field <p from <p = one 
has to invest the work 

A = -^ [ R jVv\\ 2 dV. (1.3) 



If we add the assumption that all energy acts as active gravitational mass, 
according to E = mc 2 , and also think of the integrand in (1.3) as representing 
energy density, we might be tempted to consider the modified equation 

^= teG ('-sk^» 2 ) (L4) 

with the aim to incorporate into Newtonian gravity the following 

Principle (P). All energy acts as source for the gravitational field. 

A field equation satisfying P must be non-linear. One might wonder whether 
(1.4) gives a Newtonian model that satisfies P. If it were true that it shared this 
qualitative feature with general relativity one might profitably employ this sin- 
gle scalar equation to study certain qualitative features of general relativity in a 
mathematically simpler environment. In fact, (1.4) is often proposed in pedagogi- 
cal discussions to precisely this end. For example, in [2] [3] the authors suggest that 
some useful lessons concerning the energy-regulating power of the gravitational field 
can be learned from model theories of charged particles based on (1.4). In passing 
we remark that (1.4) can be written in a linear form by making the field- redefinition 
if) := exp(</?/2c 2 ): 

. , 2irG , .„ „. 

= — pV (1-5) 
c z 

where the boundary conditions <p(r — > oo) = translate to tp(r — > oo) = 1. In the 
following we shall for simplicity always assume p to have compact support B C R 3 . 

In section 3 we discuss what is wrong with a theory based on (1.2,4) and suggest 
a different and consistent theory in section 4. Section 5 briefly discusses some 
properties of spherically symmetric solutions to the latter and section 6 points out 
the relation to general relativity. Section 2 summarizes some facts from Newtonian 
gravity. We employ the standard summation convention for repeated indices in 
up-down positions and use the euclidean metric 5 a b to raise and lower indices. V a 
denotes the partial derivative with respect to x a . 3-component vectors are also 
written with an arrow, £, with £ • f} = ^ a rj a denoting the scalar product. 
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2. Newtonian Recollections 



To see what is wrong with (1.4) it is helpful to first give a derivation of (1.3). 
Consider a one-parameter family of diffeomorphisms s h d s such that <j s= o = id 
and -^\ s= oa s (x) = |*(af). We wish to use a s to redistribute the matter by dragging 
it along this flow. Pulling back the 3-form pdV by the inverse diffeomorphisms we 
obtain p s dV := (a~ 1 )*(pdV) and hence for the Lie-derivative of the density p along 

— * 

Sp:=^ =-V-( P a (2-1) 

ds s=o 

where here and in the following we use the variational symbol, 8, for the derivative 
at s = and call it 'the variation' of the quantity in question. 

The variation of the work done to the system is easily determined using (1.2): 
SA = - f £■ fdV = - [ pV • (pi) dV. (2.2) 

JRZ JB 

Equations (2.1,2) imply 

SA= I pdpdV. (2.3) 



JB 

This equation is independent of the field equation. If we assume the validity of 
(1.1) throughout the (adiabatic) motion we can use it to eliminate dp and write 
(2.3) solely in terms of (p. The result (1.3) then easily follows. 

From (1.1,2) it follows that the force per unit mass may be derived from a 
symmetric stress tensor, f a = — V b t a &, where 

tab = ((V o¥ >)(V 6¥ >) - i<U||V^|| 2 ) (2.4) 
SA=-f U a dV=f CV b t ab dV=f V^eKbdV. (2.5) 



so that 



Here we assumed |||*(r — > oo)|| < ar for some real constant a and that derivatives 
of ip fall off as fast as r~ 2 . Vector fields which satisfy V( a £ 6 ) = (Killing equation) 
generate rigid motions and are given by f*(af) = k (translations) and £(x) = k x x 
(rotations), for constant k. For those 5 A = 0, as it must be by the principle of 
action = reaction. Otherwise the system would self-accelerate. 
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Finally we define the gravitational mass as the total flux of the gravitational 
field <p out to infinity: 



M g := lim — *— / n-Vipdo. (2.6) 
' 47rG J si 



S 2 denotes a two-sphere of radius r, n its (outward pointing) normal, and do the 
surface element on S 2 . The limit of integrals in (1.14) is sometimes abbreviated by 

is? • 



Why Inconsistent? 

Since in Newtonian theory M g = M m := f pdV, M g only depends on the amount 
but not on the distribution of matter and clearly P cannot be satisfied. Now, 
replacing (1.1) by (1.4), one obtains the following formula for the variation 5M g 

where we have used (1.4) N times to replace A(p. For a regular matter distribution 
(p will be bounded, say (p(x) < K , Va? G -R 3 . Also, the integral over j^8(A(p) 
just represents 5M g so that the last term on the right hand side is majorized by 
jjt(K I 'c 2 ) N 5M g . It vanishes in the limit N — > oo. In this limit the sum on the 
right side is just the exponential function. Thus we obtain the result: 

SM g = [ Sp exp(</Vc 2 ) dV (3.2) 
Jb 

which, recalling (2.3), deviates from 5A/c 2 by all the higher-than-linear terms in the 
expansion of the exponential. Hence (1.2,4) violates P. This is not really surprising, 
since (1.3) was derived under the assumption of (1.1,2). Changing it to (1.2,4) also 
invalidates (1.3). A correct procedure must iterate the step that led from (1.1) 
to (1.4). For example, the next (second) step would be to determine a modified 
expression for the gravitational field energy from (1.2,4) and then change (1.4) 
accordingly. Eventually this procedure should converge to a self-consistent field 
equation. However, as we will see in the next section, such a self-consistent field 
equation can actually be guessed directly. 
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At the end of this section we also point out another flaw in the combination 
(1.2,4). Using (1.5) to replace p in (1.2) one easily derives 



f a = - exp(-^/c 2 ) V & (exp(vVc 2 ) *«&) (3-3) 

with tab from (2.4). From this expression it follows that the force density is not the 
divergence of a stress tensor. There are many ways to isolate the part that obstructs 
the right hand side of (2.11) to be written as the divergence of a symmetric tensor. 
Two obvious ways are 

fa = - V b t ab - ^±^\\Vp\\ 2 V a p (3.4a) 

= - V b [(1 + p/c 2 )t ab ] + ^2 Av?V aV ? 2 . (3.46) 

The system (1.2,4) thus potentially violates the principle action = reaction. 1 



4. A Consistent Modification 

Equation (3.2) was derived assuming (1.4) but not (1.2). If we maintain (1.4,5) 
but call <p = c 2 exp(</?/c 2 ) rather than ip the gravitational potential, we have (3.2) 
just expressing the validity of P, i.e. c 2 5M g = S A with S A given by (2.3). This 
re-interpretation implies that (1.2) has to be replaced by 

/ = -PV0 (4.1) 

and that (1.4,5) written in terms of <p reads 

To be sure, for explicit calculations one would preferably use (1.5) where ip = 
c\f§. 4> must satisfy the boundary conditions 4>{r — > oo) = c 2 . The Newtonian 
approximation is obtained from expanding <p = c 2 + ip + 0(ip 2 ) and keeping only 
linear terms in (4.2). Note that in the expression (2.6) for M g we must write <fi 

1 To manifestly show a violation one should prove existence of a regular solution to (1.4) with 
tp(r — > oo) = for which f f a £, a ^ for some generator £ of a rigid motion. 
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instead of <p. But for r — > oo only the linear term in <p contributes to the surface 
integral so that (3.2) is still valid. This is why (3.2) indeed expresses the validity 
of P for (4.1,2). To be sure, once (4.1,2) are established, the equation c 2 5M g = 6 A 
is most easily proven directly. For completeness we give a short direct proof in the 
appendix. The point of our derivation of (3.2) was that it suggested the definition 
of (p in terms of (p and hence (4.2). It is interesting to note that (4.2) is precisely 
the equation that Einstein already proposed before the advent of general relativity 
in 1912 [4]. 

Equations (4.1,2) also manifestly implie the principle action = reaction in the 
sense above. Indeed, we now have f a = — /9exp(</?/c 2 )V a </?. Replacing —pV a <P by the 
right hand side of (3.3) just cancels the exponential function outside the derivative 
so that the remaining divergence can be rewritten in terms of <fi. This leads to the 
desired formula, f a = — V b t a &, with 

tab = [(V a 0)(V 6 0) - !<U||V0|| 2 ] } . (4.3) 



We may interpret the two terms on the right hand side of (4.2) as energy 
densities due to ponderable matter and the gravitational field respectively. The sum 
of both determines the convergence A<p of the gravitational field — V0. Both terms 
are positive since <p is positive. This is in contrast to (1.4), where the Newtonian 
gravitational field energy was negative definite, which is usually said to have its 
origin in the attractivity of gravity. But of course here gravity is also attractive. 
What is different here is that the (rest-) energy of matter depends on the value of 
the gravitational potential at its location. This allows that a contraction of a matter 
distribution enhances the field energy although the total energy decreases. This is 
achieved by displacing the matter into regions of smaller gravitational potential and 
thereby sufficiently decreasing the matter part of the energy. 

The total gravitational energy is given by 

E total :=c 2 M g = f p( f>dV+-^- f ^^dV =:E matter + E^ A (4.4) 

JB oTTtr J R 3 (p 

where the expression for i^id can also be written in terms of an integral over B 
(the support of p) only. To see this we recall that for large distances from the source 
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we have the expansions for cf> and 



= 1 - 



GM, 



F + 0(r" 2 ) 



err 



GM Q 9 , 



(4.5a) 
(4.56) 



so that -Etotai can also be expressed as an integral of ^gAip = c 2 ptfj (using (1.5)) 
over B. Replacing -Etotai by this expression in -Efieid = -Etotai — -^matter one obtains 



-Efield 




1 - 




dV. 



(4.6) 



5. Solution for Homogeneous Spherical Star 

In this section we determine the gravitational field for the externally prescribed 
mass distribution 



for r < R 



3M m 
47ri? 3 

(5.1) 

for r > R 

where M m is the total ('bare') mass of matter: M m = J B pdV. It will be convenient 
to introduce the 'matter radius' lZ m and the 'gravitational radius' 1Z g : 

_ GM m GMg 

C z C z 

and the abbreviation 



- = I (5-3) 

V 2R R v ' 

We use (1.5), set ip(r) = x{ r )/ r i an d obtain 

y" = J ^ X for r < R , - 4) 

A \0 for r>i?. V ' ; 

The general solution which makes (and hence ip) finite at r = is 

«->=<f-£ !:::>! (5 ' 5) 



The integration constants K and 1Z 9 are determined by the requirement that 
hence tp) should be continuously differentiable at r = R: 



TZ„ = 2R 



K 



1 - 
1 



ujR 



ui cosh(ujR) 



(and 

(5.6) 
(5.7) 



Fixing the radius R in (4.6) gives us 1Z 9 as function of lZ m , i.e., the gravitational 
mass as function of the bare mass. In terms of the dimensionless quantities y = 
TZg/R and x = 1Z m /R this reads: 



y = fix) = 2 



l 



tanh(3x/2)3 
(3x/2)* 



(5.8) 



which for x > maps monotonically [0, oo] — > [0,2]. For small x one has f(x) = 
x - |x 2 + j^x 3 + 0(x 4 ). The fact that f(x) < 2 Vx G R+ means that the 
gravitational mass is bounded by a quantity depending only on the geometry (here 
R) of the mass distribution: 

,2c 2 



M g < R-. 



(5.9) 



Note that this is achieved with all contributions to the gravitational mass on the 
right hand side of (3.3) being positive. No subtractions are taking place. Rather, 
high matter densities p are suppressed by the small potentials (p produced by them 
(i.e. 'red-shifted' in general relativistic terminology). This can be seen in detail 
from the following expressions: 



-Etotal 



2Rc 4 
G 

M m c 2 



1 - 



tanh(a;-R) 



1 - 



uR 
37?. m 



E 



matter 



Rc 4 
~G 



M m c z 



5R 
tanh(c<;-R) 
ujR 

61Zrr, 



-Afield = 



3Rc 4 



G 



M m c 2 



1 - 

1 - 
'3TZ 



5R 
tanh(a;-R) 



+ o(n 2 m /R 2 ) 

+ tanh 2 (ujR) - 1 
+ 0(K 2 JR 2 ) 



5R 



cvR 

+ o(n 2 m /R 2 ) 



- ^tanh 2 ^^) 



(5.10a) 
(5.106) 
(5.11a) 
(5.116) 
(5.12a) 
(5.126) 
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where the second expressions on the right hand sides are expansions of the first in 
terms of lZ m /R. Also, recall that 7£ m c 4 /G = M m c 2 . Note the familiar |-term in 
(5.10b) for the Newtonian binding energy. 

Decreasing R for fixed lZ m we see from (5.10b-12b) that to first approximation 
this enhances the field energy and at the same time decreases the matter energy 
twice as fast, so as to decrease the total energy by the same amount by which the 
field energy increased. Clearly the total energy must decrease in accordance with 
the attractivity of the gravitational interaction. 

Coming back to (5.9) we next show that it remains valid for any spherically 
symmetric matter distribution. In particular, it remains valid for more realistic 
matter distributions (of compact support r < R) which are determined by a coupled 
system of (4.2) with some equations of state for the matter. The proof is simply 
this: For r > R (5.4) is solved by x+( r ) = r — 1Z g /2 and by some function X-( r ) 
for r < R. Continuity and differentiability of </> at r = R is equivalent to 

x-(R) = R-\n g (5.13) 

X '-(R) = l. (5.14) 

Suppose x(-R) < 0, then x" = ^f-px with p > implies x"( r ) < for all r < R, 
with strict inequality if r lies in the support of p. Equations (5.13, 14) now imply 
that the curve r — > X-( r ) nes below the curve r — > r — ^lZ g for r < R, which in turn 
implies %(r = 0) < — \Tl g < 0, where the last inequality just expresses the positivity 
of the gravitational mass. But this contradicts the regularity of the gravitational 
potential which requires a finite value of ip(r = 0) and thus x(r = 0) = 0. Hence 
we must have x(R) > or, by (5.13), TZ g < 2R. 

Finally we mention that the finite (negative) bare mass for the point-charge 
model of ref. [V] crucially depends on taking (p (and (1.4)) rather than as gravita- 
tional potential. In the latter case it unfortunately turns out to be infinite. At this 
state of affairs the finite bare mass obtained in [V] is not a consequence of P, but 
rather an artifact of (1.4), which violates P. But P seems crucial for any model of 
general relativity and hence (1.4) does not seem suited to model general relativity 
for questions concerning the energy regulating power of the gravitational field. 
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6. Connection with General Relativity 



Let us consider a Lorentz metric in which only the time-time component of the 
metric differs from its Minkowski value. We write 

ds 2 = -2cpdt 2 + dx-dx (6.1) 

and require that for large spatial distances this approaches the Minkowski metric: 

lim </>(*) = ±c 2 . (6-2) 

Using the parameter t, the equations for a timelike geodesic curve boil down to 
(x = dx/dt) 

d?= -V0 + _1 (f • V0)f (6.3) 

which, neglecting terms oc iy/c) 2 for the moment, are just the Newtonian equations 
of motion for a point mass in the external potential <p. Setting for the moment 
20 = t(j 2 , the components of the Ricci tensor are most easily calculated with respect 
to the orthonormal tetrad: e t - = -i/> _1 <9t, e„ = d a for a = 1,2,3 (the hat over the 
indices signifies the orthonormality). We obtain: 

Rac = -^~V,«,c (6.4) 

R B = ^Aj, (6.5) 

where A is just the ordinary Laplacian d a d a and d a ip = ip :a . The scalar curvature 
then follows (summation over a): 

R = Raa~ Ru = -2^ _1 A^. (6.6) 



Let us now consider Einstein's equations 

Riiu = —g- (T^u ~ \9itvTx) (6.7) 

with an energy momentum tensor = pu^u u for stationary pressureless dust of 
rest-mass-density p and 4- velocity u = cip~ 1 d t . If we varied the Einstein-Hilbert 
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action only within the class of metrics of the form (6.1) we would obtain only the 
tt-component of (6.7), which explicitly reads: Aip = ^p-ptp. In terms of <j) we have 

This is almost (4.3) except for an additional factor of two in the G-dependence. 
But note that the boundary condition (6.2) differs from (4.5a) by a factor ^which 
implies that (6.8) and (4.2) have the same Newtonian limit. A solution of (4.2), like 
(5.5), can be easily turned into a solution to (5.6) if we multiply it by \ and replace 
G by 2G. Thus we conclude that (3.3) is essentially the time-time-component of 
Einstein's equations. Note however that we cannot solve the full set of Einstein's 
equations with the ansatz (5.1). In fact, adding the trace of the spatial part to 
the tt-part, the left hand side is zero according to (5.4, 5) whereas the right side 
is easily seen to be proportional to T^. So there is no solution, at least if the 
matter satisfies the dominant energy condition \T^\ > \T^\. Note that the trace 
of Einstein's equations (5.7), R = — ^^T, is already in contradiction to their tt- 
component, since T = T&a — Tit = = ~pc 2 so that with (5.6) the trace part 
reads Aip = —^f-ip which differs in sign from the equation above. This shows how 
Einstein's scalar equation of 1912 [4] is related to the time-time-component of the 
general relativistic tensor equation. 



Appendix 

In this appendix we wish to give a short direct proof that (4.2) satisfies P, i.e., that 
5 A = c 2 5M g . 

Using the generally valid equation (2.3), now with <j) replacing ip, we must 
eliminate p via (3.3). This is most easily done if we set cf> = c 2 ifj 2 and use (1.5). We 
obtain: 



5A = 



2tvG 



2tvG 



^5 



Aip 



dV 



n 



2nG j R3 
do 



[V>A(<fy) - (AVO <ty] dV (AA) 
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(A.2) 



Now, the conditions for large r imply that falls off as fast as 1/r 2 and Sip as fast 
as 1/r. Hence the second term in the last bracket does not contribute. Therefore 
we may reverse its sign and obtain 



8 A = 




5 [ (n-Vip)ipdV 
Js 2 




5 [ n-V(f)dV = c 2 5M g 
Js 2 
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